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ABSTRACT 

The recent detection of resolved radio emission from AG Dra by MERLIN reported 
by Ogley et al. is discussed in the context of the wind environment and the physical 
parameters and geometry of this symbiotic binary system. In particular, it is shown 
that the two radio components are closely aligned with the binary axis, and their 
separation suggests their origin in jets ejected from AG Dra during the recent 1995- 
98 series of oubtbursts. 
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1 INTRODUCTION 

AG Dra is the well-studied symbiotic binary consisting of 
a high-velocity, metal-poor, bright K-type giant and a hot 
white dwarf companion (e.g. Mikolajewska et al. 1995; Smith 

\ et al. 1996). The binary system has an orbital period of 549 
days and a well-defined spectroscopic orbit (Fekel et al. 2001, 
and references therein). AG Dra is also among the most ac- 

, tive symbiotic stars. Its optical light curve is characterized 
by a series of active (outbursts) and quiescent phases (e.g. 
Fig. 7 of Galis et al. 1999). Although the activity of AG 
Dra, and other classical symbiotic stars is still poorly un- 
derstood, multifrequency observations covering a few whole 
activity cycles indicate that in AG Dra this activity is related 
to changes of both radius and temperature of the hot com- 
ponent (e.g. Mikolajewska et al. 1995; Greiner et al. 1997). 
AG Dra is also one of 10 known galactic supersoft X-ray 
sources (Greiner et al. 1997). 

Recently, the star has been searched, together with all 
the northern supersoft X-ray sources, for radio emission at 
5 and 8.4 GHz (Ogley et al. 2002, hereafer O02). The ob- 
servations by MERLIN telescope have confirmed previous 
VLA detections. Moreover, the source has been resolved at 
the milliarcsec scale into two components of nearly equal 
brightness, and combined flux of ~ 1 mjy. O02 have also 
studied possible interpretations of this emission in terms of 
a wind environment from either the cool giant or the hot 
white dwarf. They have concluded that all their scenarios 
give the radio emission fluxes an order of amplitude lower 
than the observed value. 

In the following we reanalyze a possible origin of the 
resolved radio emission from various wind environment, and 
show that it presumably arises from jet(s) ejected from the 
hot component during its recent series of outbursts. 



2 RADIO EMISSION FROM AG DRA 

2.1 Spherically symmetric steady outflow 

We point out here that O02 have adopted unreasonably low 
values for both the cool giant wind and the hot component 
luminosity in their analysis. Below in this section, we provide 
a critical analysis of their assumptions and calculations. 

The wind environment of AG Dra has been discussed 
in several papers, and estimates for the mass loss from 
both components are available. All these estimates have 
been based on the assumption of an isothermal, spheri- 
cally symmetric steady flow. In particular, the mm-submm 
radio spectrum and K — [12] colour excess are both con- 
sistent with the cool giant mass- loss rate, Mg/v ~ 2.5 x 
10~ 8 (d/2.5kpc) 3/2 M Q yr~ 1 /( kms " 1 ) (Mikolajewska, Ivi- 
son & Omont 2002; Kenyon 1988), whereas the UV con- 
tinuum and emission line analysis (Miirset et al. 1991) gave 
M g /w>4 x 10~ 8 (d/2.5kpc) M Q yr- 1 /(kms- 1 ), in all cases 
values two orders of magnitude larger than M g /i> ~ 5 x 
10~ 10 M Q yr _1 /( kms_1 ) adopted by O02. Using the Wright 
& Barlow (1975) formula for completely ionized wind, the 
first estimate above gives a flux at 5 GHz of ~ 0.5 mjy inde- 
pendent of the adopted distance. Thus, the cool component 
wind could, in principle, account for the observed intensity 
of the radio emission but not for the double structure. 

We also find an error in the value of irradiation luminos- 
ity adopted by O02 in their estimate for the mass loss rate 
in the cool giant due to irradiation by its hot companion, 
L h = 1.4 x 10 36 ergs -1 (or L 30 = 0.14 in units of 10 30 Js" 1 ; 
their Eq. (3) and Table 2). This error apparently stems from 
a mistake in rescaling results of Greiner (2000) to d = 1.7 kpc 
used by O02. The rescaling gives in fact L30 = 0.44. Using 
Eq.(3) of O02 with more realistic values, namely L30 ~ 1 
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Table 1 . History of radio observations of AG Dra 



JD 


Frequency 


Flux 


Reference 


2400000+ 


[GHz] 


[mjy] 




45006 


4.9 


< 0.41 


ST90 


46515 


4.9 


0.60 ±0.21 


ST90 


46611 


4.9 


>0.5 


TC87 


46646 


4.9 


0.36± 0.08 


ST90 


46646 


14.9 


0.77 ±0.23 


ST90 


48290 


8.3 


< 0.17 


SKT93 


51622-30 


5.0 


~ 1 


O02 



ST90 - Seaquist & Taylor 1990; SKT93 - Seaquist et al. 1993; 
TC87 - Torbett & Campbell 1987. 

(e.g. Greiner et al. 1997), r 2 ~ 80 Rq, and d ~ 2.5 kpcj] 
(Mikolajewska et al. 1995), we find the mass loss rate for the 
evaporated wind, M g ~ 1.2 x 10 -7 Mq yr _1 , and the pre- 
dicted quiescent 5 GHz flux of ~ 0.1 (d/2.5 kpc) -2 / 3 mjy. 
This radio flux can be further increased by increasing the 
hot component luminosity (e.g. Mikolajewska et al. (1995) 
found an order of magnitude increase in Lh due to activity) , 
however, this model still does not account for the double 
source structure. 

Concluding this section, the resolved radio emission 
from AG Dra detected by MERLIN cannot be interpreted 
in terms of the spherically symmetric steady wind from the 
cool giant. We note, however, that the Hipparcos position of 
AG Dra (Perryman et al. 1997) practically coincides with 
that of the Nl component in the MERLIN image. This in- 
dicates that the Nl component can, in principle, originate 
in the cool giant wind. 

2.2 Episodic flow 

The active phases of AG Dra are characterized by an in- 
crease of the emission line widths and marked P Cyg struc- 
ture of high ionization UV lines (e.g. Nv), which indicate 
outflow velocities of 200-300 km s _1 (Mikolajewska et al. 
1995, and references therein), while the optical Hell and 
Hi Balmer lines develop broad emission wings (e.g. Tomova 
& Tomov 1999). This behaviour suggests that the hot com- 
ponent develops a significant wind in outburst. Tomova & 
Tomov (1999) estimated the wind velocity of ~ 800kms _1 
and the mass loss rate of ~ 2 x 10 -7 Mq yr _1 from the broad 
wings of Ha and H/3 profiles observed during the 1995 out- 
burst. Their estimate made use of the distance-dependent 
estimate for the hot component radius of Mikolajewska et 
al. (1995), and thus it corresponds to d = 2.5 kpc. It was 
also based on the assumption of the spherically symmetric 

1 Parenthetically, we also point out that the values of d = 1.7 kpc 
and T2 = 30 Rq adopted by O02 from Tomov, Tomova & Ivanova 
(2000) are mutually inconsistent given the observed near IR mag- 
nitudes, colours and spectral classification of the giant. Namely, 
if we accept r 2 = 30 Rq, we obtain (using the Barnes-Evans re- 
lation, Cahn 1980, with the surface brightness appropriate for a 
mid K-typc giant and the observed K = 6.2, e.g. Belczyhski et al. 
2000) d ?s 1 kpc, even more discrepant with our preferred value 
of d = 2.5 kpc. 



steady wind which is probably not true for the hot compo- 
nent wind developed during the outburst. Anyway, applying 
the Wright & Barlow formula (1975), such a wind should 
give rise to the radio flux of ~ 5 /iJy at 5 GHz, again much 
lower than the flux observed by O02. 

Although only a few attempts have been made to ob- 
serve AG Dra at radio wavelenghts, it seems that the radio 
emission may be variable, and possibly related to the hot 
component activity (Table 1). AG Dra was not detected at 
4.9 GHz on JD 2 445 006 (Seaquist & Taylor 1990) when 
the star was at maximum of a large eruption (see the light 
curves in Mikolajewska et al. 1995, and Galis et al. 1999), 
whereas it was detected at 3 epochs (Torbett & Campbell 
1987; Seaguist & Taylor 1990) during less pronounced burst 
and its decline a few years later. It is particularly important 
that AG Dra was not detected on JD 2 448 290, when the 
star was in quiescent phase, by the highly sensitive VLA 
survey at 8.3 GHz (Seaquist et al. 1993), and the 3<r upper 
limit for the 8.3 GHz flux of <0.17mJy was much lower 
than the 5 GHz fluxes of any of the two radio components 
detected by O02. 

The hot component entered another large outburst in 
1995, which was followed by a series of more or less promi- 
nent bursts of activity till at least the end of 1999. It is 
interesting that there are two kinds of outbursts: stronger, 
cool outbursts during which the hot component tempera- 
ture decreases, and fainter, hot outbursts during which it 
increases (Gonzalez-Riestra et al. 1999; Mikolajewska et al. 
1995). The temperature decrease during the strong cool out- 
bursts is very likely due to slow expansion of the hot com- 
ponent to about 2-3 times its original size (Mikolajewska et 
al. 1995; Greiner et al. 1997) and possible development of an 
optically thick wind. All the previous radio detections were 
made during hot outbursts. Moreover, in 1997 and 1998, 
AG Dra was detected at 1.3 and 0.85 mm by the IRAM and 
JMCT telescopes, respectively (Mikolajewska et al. 2002). 
At the time of both detections AG Dra was at the hot ac- 
tivity stage. 

The fact that the radio emission is resolved into two 
compact components raises an interesting possibility that it 
originates in jets ejected from AG Dra. Bipolar jets have 
been detected in a few symbiotic systems, e.g. in V694Mon, 
CH Cyg, Hen 3-1341 and StHA 190 (Belczynski et al. 2000, 
and references therein; Tomov, Munari & Maresse 2000; Mu- 
nari et al. 2001), and in all cases they have been associated 
with the hot component outbursts and an appearance of 
a high- velocity stellar wind. A relation between the radio 
emission and the hot component activity stage has been also 
found in another symbiotic binary, CI Cyg. The cm-submm 
spectrum of CI Cyg is inconsistent with predictions of any 
model involving spherically symmetric wind(s), and the ra- 
dio emission most likely originates from a bipolar outflow 
(Mikolajewska & Ivison 2001). 

Figure 1 shows the relationship between the binary ge- 
ometry and the resolved MERLIN image of AG Dra. From 
the data in Table 1 of O02, we have estimated a separa- 
tion of the components, a ~ 0.4 ±0.1 arcsec, and a po- 
sition angle of the whole structure, PA ra d « 35 ± 20°. 
From spectropolarimetry we know the binary orbit orien- 
tation, Q — 150 ± 20°, as well as the orbit inclination, 
i = 120 ± 20° (Schmid & Schild 1997). Thus the two ra- 
dio components are within the observational errors aligned 
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Figure 1. Schematic illustration of the relationship between the 
binary orbit, derived from spcctropolarimetric observations of the 
A 6825 feature by Schmid & Schild (1997), and the Nl and N2 
features of the resolved MERLIN image of AG Dra at 4.994 GHz 
in 2000 March (O02). The orbit of the hot component relative to 
the cool giant (marked as the grey disk) is given by dots (not to 
be scaled relative to the MERLIN image). The orientation of the 
orbital plane is S7 fs 150° , and the system inclination is i 120° . 



with the binary axis (fi - PA rad ~ 115 ± 30°; Fig. 1). 
Moreover, the possible extended emission reported by Tor- 
bett & Campbell (1987) has similar orientation. The sepa- 
ration of the two radio components resolved by MERLIN 
corresponds to As ~ 1000 [d/2.5 kpc] a.u. Assuming that 
they have been ejected perpendicularly to the orbital plane 
with a velocity of the order of the escape velocity from 
the hot component, v c = (2GM h /R h ) 1/2 ~ 800 km s" 1 for 
M b ~ 0.4-0.6 Mq, and R b ~ 0.22 Rq at the maximum of 
the large otbursts (Mikolajewska et al. 1995; Greiner et al. 
1997), this separation implies that the ejection took place 
At ~ As/2w e cosi ~ 3 years before the MERLIN observa- 
tions, possibly during the recent 1995-98 series of outbursts. 



3 CONCLUSIONS 

The major results and conclusions of this paper can be sum- 
marised as follows: 

(i) The radio emission from the symbiotic binary AG Dra 
seems to be variable, and probably related to the hot com- 
ponent activity. 

(ii) The cool component wind can, in principle, account 
for the intensity of the Nl component (O02), which position 
practically coincides with the Hipparcos position of AG Dra. 

(iii) The two radio components resolved by MERLIN 
(O02) are practically aligned with the binary axis of AG Dra. 
The possible extended radio emission reported by Torbett 
& Campbell (1987) has similar orientation. The resolved ra- 



dio emission presumably originates in jets ejected from the 
binary system. 

(iv) Assuming that the jet velocity is of order of the escape 
velocity of the hot component, the separation between the 
two radio sources indicates that the ejection took place ~ 3 
yr earlier, and it was associated with the recent series of 
outbursts. 
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